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ABSTRACT 

A topology model constructed from surface-streamer visualisation describes the flow around a 
generic conventional submarine hull form at pure yaw angles of 0°, 10° and 18°. The model is 
used to develop equations for sway-force and yaw-moment coefficients which relate to the hull- 
form geometry and the flow circulation around the submarine. 
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Topology Model of the Flow around 
a Submarine Hull Form 


Executive Summary 

An understanding of the flow around a submarine is necessary to predict the hydrodynamic forces 
which can affect motion control, safe manoeuvring and signature performance of the submarine. 

The focus of this study is on the flow produced by a submarine at yaw angles of 0°, 10° and 18 °. 
Notable flow structures are the counter-rotating longitudinal vortices leeward of the hull and the 
tip vortex aft of the fin (or sail). 

In this report, a topology model of the flow around the submarine has been constructed by using 
flow-visualisation data obtained from (i) testing with a generic conventional “Joubert” hull form 
in the Defence Science and Technology Group low-speed wind tunnel at Fishermans Bend, and 
from (ii) literature review of classical structures known to exist in submarine flows. 

The topology model provides an initial step towards developing equations which relate the sway- 
force and yaw-moment coefficients to the submarine geometry and its surrounding flow. The 
model equations take into account the slenderness ratio of the hull and the location of the fin, and 
therefore may be used to compare different generic submarine designs. The equations are cali¬ 
brated by using circulation, force and moment data from experiments and/or computational-fluid- 
dynamics (CFD) modelling, and so it is possible to use them for cross-checking measurements 
and for validating CFD results. 
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Flow-topology (bifurcation and critical point) symbols 

B + positive-bifurcation line 

B negative-bifurcation line 

F foe us/vortex 

N node 

S saddle point 

Geometric parameters and body coordinate system 

A/,(x) cross-section area of the submarine hull at a given body-axis (x) location 
A p planform area of the submarine bare hull 
Cf chord length of the tin 

h c height of the casing measured from the body (x) axis 
hf height of the fin measured from the casing 

L length of the submarine 

L c /4 distance from the submarine nose to the quarter-chord point of the (NACA) fin 

Lf distance from the submarine nose to the leading edge of the tin 

L m distance from the tail plane to the mid-ship (or center of buoyancy) of the submarine 

R slenderness ratio of the submarine bare hull, R = L/(!r m ) 

r m radius of the bare hull at mid-ship (or maximum radius of the bare hull) 

x, y, z body-axis, horizontal and vertical directions 

Flow parameters 

Cf sway-force coefficient 

Cm. yaw-moment coefficient 

F y j sway force acting on the fin 

F y h sway force acting on the hull form 

i y ji(x) hydrodynamic impulse in the cross-stream (y) direction of the hull form 

K\, K 2 non-dimensionalised moments of vorticity 

Kf, Kj, model coefficients for the load on the fin and on the hull form 
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M z j\ m yaw moment acting on the tin about the mid-ship of the submarine 
yi Zt h,m yaw moment acting on the hull form about the mid-ship of the submarine 
Rei Reynolds number based on the length of the submarine, Re/ = LLL/v 
Ho free-stream velocity 

z c vertical location of the centroid of a vortex (or region of concentrated vorticity w A ) 

Greek symbols 

r circulation 

y f model coefficient for vertical location of the fin-tip vortex 

y/j model coefficient for vertical distance between the counter-rotating hull vortices 

K dimensionless circulation parameter, K = T/ [r m LL sin(\|/)] 

V kinematic viscosity of the working fluid 

CO* body-axis (x) component of vorticity 

\|t submarine yaw angle 

p density of the working fluid 

Subscript symbols 

a aft-body (or stern) of the submarine 

c casing (top side) 

fb fin bound vortex (leeward side, for yaw conditions) 
fj fin junction (leading edge) 

fl fin leeward side (for yaw conditions) 

fp, fs fin port and starboard (for straight-ahead and yaw conditions) 
ft fin tip (trailing edge, for yaw conditions) 
hb, he hull base and casing (leeward side, for yaw conditions) 
hi, hw hull leeward and windward sides (for yaw conditions) 
n nose (or bow) of the submarine 
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1 Introduction 


A submarine performing a yaw manoeuvre can produce some noticeable large-scale flow struc¬ 
tures, such as a pair of counter-rotating longitudinal vortices leeward of the hull and a tip vortex aft 
of the tin. The circulation of these structures are of concern because they can affect the hydrody¬ 
namic forces acting on the hull, and hence the motion control and computational-fluid-dynamics 
(CFD) prediction of safe manoeuvre. Also another major factor is the vibration and noise due 
to interactions between flow structures and the submarine propulsion, which can affect signature 
performance. 

Low-speed wind-tunnel testing and CFD modelling of submarines by the Defence Science and 
Technology Group at Fishermans Bend make use of a generic conventional hull form 1 defined by 
Joubert [1,2]. The design objectives of this hull form are to provide “ mi nimum practical resistance 
and minimum water flow noise... while still carrying out all its normal functions” [2]. To reduce 
resistance and flow noise arising from flow-structure interaction, it is necessary to test the shape 
of the submarine, which includes the length-to-diameter (i.e. slenderness) ratio of the hull and 
the location and height of the fin 2 . A key aspect of design testing and CFD validation therefore 
requires an understanding of the flow structures produced by different parts of the submarine and 
how they contribute to the hydrodynamic forces during a yaw manoeuvre. 

The study presented in this report is twofold: (i) to visualise and establish a detailed topology of 
flow structures produced by different parts of the submarine (Sections 3-7), and (ii) to develop 
equations which relate the hydrodynamic forces to the circulation of the flow structures and the 
shape of the submarine (Sections 8-10). 


2 Shape of the Submarine Model 


Figure 1 shows the shape of the generic hull form used for this study. The model is developed 
from a bare hull which has a slenderness ratio of 7.3 to provide minimum resistance [2], The nose 
geometry is derived from a NACA-0014.2-N00.20 profile and is axisymmetric for the first 7% 
of the body length (L) to provide a clean forward sonar. The fin has the shape of a NACA-0015 
aerofoil 3 with a rounded trailing edge; the fin height is 8%L, the chord length is 16 %L and the 
leading edge is located at 31%L. The tapering to the end of the tail cone starts at approximately 
76%L. The “X”-rudders located at approximately 86-91 %L are the aft control surfaces. The 
present model does not include any hydroplanes. 


'The term “hull form” used here takes into account the profile of the submarine nose (or bow), the casing and the 
aft-body (or stem); an axisymmetric hull form with no casing is called a bare hull. 

2 The fin (Commonwealth and European) or the sail (American) refers to the largest appendage on the upper hull of 
the submarine geometry. 

3 The NACA series is a well-known family of aerofoil shapes developed by the American National Advisory Com¬ 
mittee for Aeronautics — NACA has published a wealth of information on aerofoil lift and drag coefficients established 
from wind-tunnel testing, e.g. [5], 
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NACA-0015 Fin 



Figure 1: A schematic diagram of the generic conventional hull form defined by Joubert [2 ]. 


3 Surface-Streamer Visualisation 

Evidence of strong turbulence fluctuation or flow separation is obtained by attaching a series of 
streamers on the submarine model of length L = 1.35 m. The streamers are woollen tufts approx¬ 
imately 10 mm long and 1 mm in diameter. Motion of the streamers in the air flow is observed in 
the test section (2.7 m wide x 2.1 m high x 6.6 m long) of the low-speed wind tunnel at Fisher¬ 
mans Bend [6]. To ensure that the surface shear stress is high enough to make the tufts flutter and 
that the results are compatible with earlier experiments and simulations [7, 8, 9], the test is per¬ 
formed at a free-stream velocity (LL) of 50 m/s. The Reynolds number based on the body length 
of the submarine is sufficiently large (Rc[ = LLL/v = 4.5 x 10 6 ) such that details of the flow are 
not sensitive to small changes in the Reynolds number. 
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To adequately capture details of the surface flow, the streamers are closely distributed on the 
submarine; to reduce setup time, they are attached only on one side of the body. With the subma¬ 
rine supported on a turntable via two streamlined pylons, it is possible to yaw the body; surface- 
flow patterns on the windward and leeward sides are obtained by yawing towards starboard and 
port, respectively. Top-view and side-view images are obtained for three different yaw angles 
\jt = 0°, 10° and 18°. Each image is an ensemble average of 125 frames from a digital video 
camera (25 frames per second) with the background subtracted to isolate the streamers. 


4 Topology of the Submarine Flow 

For straight-ahead condition (at zero yaw), the surface-flow pattern is symmetrical about the mirror 
plane of the submarine (Fig. 2). Inspection of video recording has indicated that the flow pattern is 
stable, and that there is no visible large-scale separation along the hull. The simplest interpretation 
in Fig. 3 shows that the surface streaklines begin at a reattachment node (N n ) on the nose and 
terminate at a separation node (N a ) aft of the submarine. Notable junction-flow features produced 
by the fin on the casing include a stagnation node-and-saddle (Nfj and Sfj) pair and a “U-shaped” 
negative-bifurcation (B~) line. On the casing, a positive-bifurcation (B ( l ) line runs between the 
legs of the U-shaped negative bifurcation ( B ~). The enclosed surface-flow pattern (Fig. 3) satisfies 
the spherical topology [3, 4]: 

number of nodes — number of saddles = 2. (1) 


At 10° yaw, the flow over the submarine is no longer symmetrical (Fig. 4). As the surface of the 
hull is continuous in the circumferential direction, Fig. 5 shows that the flow spreading from both 
the windward side (positive bifurcation B7) and the leeward side (positive bifurcation B h! ) con¬ 
verges on the upper hull (negative bifurcations BJ u , and B~) and on the lower hull (negative bifurca¬ 
tion Bj ib ), thus satisfying continuity. On the leeward side of the fin, a streak of fluctuating streamers 
indicates strong turbulence and separation. In Fig. 5, this feature is shown as a negative-bifurcation 
(Bji) line. Further inspection of video recording have shown possible counter-clockwise rotation 
of streamers at the leading edge of this bifurcation (B^) line; in the alternate interpretation (Fig. 5), 
this supplementary feature appears as a surface vortex (focus Fjf, the existence of saddle S*p is de¬ 
duced from Eq. (1)). Note that the flow in this region is highly turbulent, and that the surface 
vortex (Bp) may be intermittent. 

At 18 ° yaw, the surface flow is different from that obtained at 10° yaw; namely, the flow gradients 
are steeper along the lines of positive bifurcation (B bw and B hl ) and there is strong turbulence 
fluctuation over a larger leeward area of the fin (Fig. 6(a)). From instantaneous images and video 
recording, clearly observed streamers diverging from a point source on the leeward surface of the 
fin implies that there is a reattachment node here (Fig. 6(b)). In the interpreted flow field (Fig. 7), 
strong shear-layer interaction in the vicinity of this reattachment node (N/i) leads to a classical 
U-shaped separation (negative bifurcation Bp) around this node. For the 18°-yaw case to be 
consistent with the observations for the 10°-yaw case, it is possible to include a supplementary 
focus-and-saddle (Fp and S*p) pair in the vicinity of the U-shaped separation, where the topology 
still satisfies Eq. (1). 
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Nose (n) 


Aft-body (a) 



Figure 2: Time-averaged surface-flow pattern on the submarine; \]/ = 0°. This flow pattern is 
stable and has the same appearance as the instantaneous surface-flow pattern. The 
mean flow direction is from left to right. 




N n + N a + Nfj = £ nodes = 3 
Sfj = £ saddles = 1 


Figure 3: Interpretation of Fig. 2 using spherical topology £ nodes — £ saddles = 2. 
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(a) Time-averaged flow 




Fin leeward {ft) 


Leeward 

side 



Hull base (hb) 


(b) Instantaneous flow 



Figure 4: (a) Time-averaged and (b) instantaneous surface-flow pattern on the submarine (in 3rd- 
angle orthogonal projection); \J f = 10°. The mean flow direction is from left to right. 
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Figure 5: Streakline interpretation of Fig. 4 using spherical topology Enodes — E saddles = 2. In 
the alternate interpretation, additional critical points are marked with an asterisk (*). 
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(a) Time-averaged flow 





(b) Instantaneous flow 



Figure 6: (a) Time-averaged and (b) instantaneous surface-flow pattern on the submarine (in 3rd- 
angle orthogonal projection); \\f = 18°. The mean flow direction is from left to right. 
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N n +N a + N/j + N/i + (Ffr) = Enodes = 5 
S fj + Sfi + (S* fI ) = E saddles =T 


Figure 7: Streakline interpretation of Fig. 6 using spherical topology Enodes — E saddles = 2. In 
the alternate interpretation, additional critical points are marked with an asterisk (*). 
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From three-dimensional reconstruction of the streakline patterns in Figs. 3, 5 and 7, it is possible 

to infer the main features of the separated flow based on rules of topology [3, 4]. These flow 

features, as shown in Figs. 8, 9 and 10 for increasing angle of yaw, are identified as follows. 

Separation over the nose — As the angle of yaw is increased, the shear layer over the nose 
separates on the leeward side and rolls up into a pair of counter-rotating eddies (Figs. 8, 9 
and 10). Interaction between the shear layer and the hull/casing allows the leeward eddies 
to grow in size while the local free-stream stretches the eddies into the form of longitudinal 
vortex legs (/j, c and /•/,/,). The leeward vortex legs are on opposite sides of the positive- 
bifurcation (Bjjj) line and they run approximately adjacent to the negative-bifurcation lines 
along the hull (i.e. vortex leg F/ u: with line B h -B ( ; vortex leg F/,/, with line Bj ib ). Further 
details of the flow around the nose are given in Section 5. 

Stagnation at the fin junction — On the casing, the turbulent boundary layer separates at the 
leading edge of the fin and rolls up to form junction eddies. Adverse pressure gradient near 
the junction leads to the stretching of the eddies into the form of a “horseshoe vortex”. In 
Figs. 8, 9 and 10, the leading edge of the horseshoe vortex (Fj\ ) is located in a recircula¬ 
tion region between the stagnation points (Nfj and Sfj) at the fin junction. The legs of the 
horseshoe vortex ( Ff p Fj\Ff S ) are on opposite sides of the fin and they run adjacent to the 
U-shaped negative bifurcation (B“); increasing the yaw angle increases the asymmetry of 
the horseshoe vortex. For simplicity, the present topology excludes junction-flow features 
produced by the control surfaces and supporting pylons; for more details on junction flow, 
see Section 6. 

Formation of the fin-tip vortex — At 0° yaw (Fig. 8), the flow is symmetrical about the mir¬ 
ror plane of the fin and the separation occurs along the tip of the fin to produce a pair of 
counter-rotating longitudinal eddies. Increasing the angle of yaw reduces the symmetry of 
the separated flow and redistributes the developing eddies toward the leeward (suction) side. 

At 10° yaw (Fig. 9), an asymmetrical flow is established and the separation along the tip 
of the fin is accompanied by a strong negative bifurcation (B^) on the leeward face of the 
fin. With a sufficiently high shear stress, this bifurcation (B^) line can terminate as a focus 
(Fji). In Fig. 9, this focus (F'p) is the end of a line vortex which rises from the surface of the 
fin and, because the line vortex is embedded in the wake of the fin, it is drawn downstream 
of the fin as a longitudinal vortex. Since the vortex legs from the surface focus (F^) and 
from the tip of the fin have the same rotation sense, both legs combine to form the fin-tip 
vortex (Ff t ). 

By increasing the yaw angle from 10° to 18 °, the leeward flow further produces a reattach¬ 
ment node ( N/i ) on the fin. In Fig. 10, the shear layer separating around this node (Nfi) rolls 
up into a U-shaped vortex. Near the base of the fin, the U-shaped vortex leg combines with 
the horseshoe leg from the junction vortex (Fji) of the same rotation sense (clockwise Ff p ); 
near the top of the fin, the U-shaped vortex leg combines with the fin-tip vortex of the same 
rotation sense (counter-clockwise Ff t ). The portion of the U-shaped vortex in the spanwise 
direction of the fin is known as the “bound vortex”. 

As shown in Figs. 9 and 10 for the 10°-yaw and 18°-yaw cases respectively, the exclusion 
or inclusion of the surface vortex (F^) on the fin produces only minor changes to the overall 
flow topology. Further description of the fin-tip vortex is provided in Section 7. 


UNCLASSIFIED 


9 



DST-Group-TR-3177 


UNCLASSIFIED 


U-shape 
separation 
over casing 

b: 



Fin-tip 

vortices 


Fin-junction 
(horseshoe) vortices 



Figure 8: Three-dimensional representation of the flow field; »[/ = 0°. (a) Time-averaged surface 
streaklines and (b) vortex lines. The mean flow direction is from left to right. 
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Figure 9: Three-dimensional representation of the leeward flow field; \\f = 10°. (a) Time-averaged 
surface streaklines, (b) vortex lines and (c) alternate interpretation with surface vortex 
(Fp) on the fin. The mean flow direction is from left to right. 
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Figure 10: Three-dimensional representation of the leeward flow field; \|f = 18°. (a) Time- 

averaged surface streaklines, (b) vortex lines and (c) alternate interpretation with sur¬ 
face vortex (Ff) on the fin. The mean flow direction is from left to right. 
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5 Flow around the Nose 


Figure 11 shows the effect of increasing yaw angle on the surface streaklines and the cross-stream 
around the nose; three-dimensional representation of the flow is given in Fig. 12. The sequence of 
flow is inferred from present laboratory observations, previous work on similar ( ellipsoidal ) nose 
forms [4] and from standard topology rules. 

For zero and very small angles of yaw (0 < \(/ < 5 °), Figs. 11(a) and 12(a) show that the flow is 
everywhere attached to the nose and that surface streaklines originate from the tip of the nose at 
the node N n . With a minor asymmetry in circumferential pressure, streaklines tend to diverge on 
the windward side to produce a weak positive bifurcation (B^ w ). On the leeward side, the flow 
pressure is lower and the converging streaklines form a weak negative bifurcation (B^). In cross- 
section of the surface flow, the positive bifurcation (B^) is a “reattachment” half-saddle (S/ JVV ) 
and the negative bifurcation (B hj ) is a “separation” half-saddle (S' hl ). The interpreted flow on a 
two-dimensional plane cutting the submarine body satisfies the cross-sectional topology [3, 4]: 

number of nodes — number of saddles = — 1. (2) 


As the yaw angle is increased (5°<\|r< 10°), the flow over the nose becomes locally unsta¬ 
ble and this produces an array of dissipative structures. On the surface of the nose shown in 
Figs. 11(b) and 12(b), they appear as an array of very fine alternating local positive- and negative- 
bifurcation lines. These structures scale with the thickness of the boundary layer and they resemble 
an array of eddies moving approximately in the direction of the local mean flow. Since the struc¬ 
tures are very small, they have negligible effect on the global flow pattern; however, the inclusion 
of these (x-number of) eddies has to preserve topology (Eq. (2); Fig. 1 1(b)). 

For larger yaw angles (10° <\|/<20°), the flow near the tip of the nose becomes less stable and this 
leads to the formation of new structures similar to those developed further downstream at smaller 
angles of yaw. Downstream of the nose, the older (eddy) structures coalesce to produce larger 
local (vortex) structures. Eventually, the effect of these larger local structures becomes significant 
enough to redistribute the surface shear stress and alter the flow pattern in the vicinity of the nose. 
In cross-section of the altered surface-flow pattern (Fig. 11(c)), the “reattachment” half-saddle 
(Shi) corresponds to a positive bifurcation (B^) and the “separation” half-saddles (She and .S'/,/,) cor¬ 
respond to respective lines of negative bifurcation (BJ XC and B hh ). In Fig. 12(c), the flow separating 
from the lines of negative bifurcation (B, and £U) rolls up to produce a pair of counter-rotating 
longitudinal vortices (F) K . and /-/,/,). From visualisation experiments, these vortices have been iden¬ 
tified by helical twisting of dye (or smoke) filaments in the separated flow [10], as shown by the 
example in Fig. 13. 

With further increase in the angle of yaw (\|/>20°), the redistributed shear stress is sufficient to 
produce a pair of counter-rotating surface foci (F* lc and Fj* b ) at the leading edge of the negative- 
bifurcation (Bj tc and B /th ) lines. In Figs. 1 1(d) and 12(d), the vortex core rising from each focus is 
stretched downstream of the nose to form a vortex pair (F^ c Fh c and F^Fhb), known as the Werle- 
Legendre horn vortices. Figure 14(a) shows a classical side-view profile of a horn vortex obtained 
from dye visualisation [4]. In Fig. 14(b, c), the surface-flow pattern produced on the leeward side 
of the nose/hull is topologically the same as those shown in Figs. 11(d) and 12(d). 
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(a) 

£ nodes =0 

\Sh w + kS' hI = Esaddles = 1 

(b) 

•x(Feddies) = £ nodes =x 

2 ^hw 2^hi 2^eddies) — XSclddlcS— 1 + X 

(c, d) 

Flic + Fi,b = £ nodes = 2 

2 $hw + 2^hl ^hi 2^hc + 2 ^hb ~ ^ Saddles — 3 


Figure 11: Streaklines on the submarine nose (top view, the mean flow direction is downwards) 
showing the effect of increasing yaw angle, after [4], In cross-section (A-A), the 
flow is viewed from the downstream end of the hull, with cross-sectional topology 
£ nodes — Esaddles = — 1. For a three-dimensional representation of each surface flow, 
see Fig. 12. 
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(a) 
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Figure 12: Vortex lines on the leeward side of the nose showing the effect of increasing angle of 
yaw. For evidence of the nose/hull vortices and the Werle-Legendre focus of separation, 
see Figs. 13 and 14, respectively. The mean flow direction is from left to right. 
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Figure 13: Flow visualisation of the nose/hull vortices for \|/ = 15°; (a) dye filaments in water 
[10] and the interpretation of (b) surface streaklines and (c) vortex lines. From side 
view, the mean flow direction is downwards. In cross-section (A-A), the flow is viewed 
from downstream of the nose. 



Figure 14: Flow visualisation of the Werle-Legendre focus of separation for \| f>20°; (a) dye 
filaments in water [4] and the interpretation of (b) surface streaklines and (c) vortex 
lines. From side view, the mean flow direction is downwards. In cross-section (A-A), 
the flow is viewed from downstream of the nose. 
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6 Junction Flow at the Root of the Fin 

Figure 15 shows an enlarged view of the junction flow at the tin. In the vertical symmetry plane 
of the fin, the surface node Nfj and saddle Sfj are projected as half saddles S' ^ and Sj\ , respec¬ 
tively. In the flow approaching the fin junction, the boundary layer separates and develops a series 
of horseshoe vortices under a streamwise pressure gradient; these horseshoe vortices are subse¬ 
quently stretched around the fin under a spanwise pressure gradient. During vortex stretching, new 
vortices are formed upstream of the oldest horseshoe vortex; topology study (e.g. [11]) has shown 
that it is possible to have up to three primary horseshoe vortices injunction flows. In Fig. 15, the 
primary vortices (Fj\, Fp_ and Fp) have the same rotation sense as does the vorticity of the ap¬ 
proaching boundary flow; the accompanying eddies (F e \, F e 2 and Fp) are in the opposite rotation 
sense to satisfy the regional topology [3, 4]: 

number of nodes — number of saddles = 0. (3) 

At large Reynolds numbers, the primary vortices can induce velocity components on themselves 
and on the accompanying eddies, and this can lead to (i) the stretching or diffusion of vorticity, 
(ii) the “leapfrogging” of vortices between one another, and (iii) the merging of vortices to form 
a larger vortex [11]. Over time, the largest/oldest vortex (closest to the fin) breaks down (due to 
vortex stretching) while the newer vortices grow and subsequently combine with the remains of 
the oldest vortex; this process is repeated but is generally aperiodic [11], 

Figure 16 shows a classical time-averaged surface-oil flow pattern around a NACA fin which re¬ 
sembles topologically the streamer pattern obtained from the present geometry (Section 4). More¬ 
over, the oil-visualisation pattern in Fig. 16 has been obtained from a much larger tin geometry 
(30.5 cm chord x 22.9 cm span) [12], which makes it possible to observe the finer details of the 
surface flow — for example, the line of low mean shear between the primary stagnation point and 
the tin. 

In the enlarged view provided in Fig. 17, the regions of low mean shear stress and/or stagnation are 
indicated by areas with high deposit of oil (or large patches of white). From Fig. 17, it is possible 
to describe the primary horseshoe vortices in relation to the surface-flow pattern upstream of the 
fin. The interpreted streaklines in Fig. 18 show the primary stagnation as a saddle point (Sfj)\ the 
surface flow separating from this saddle point (Sfj) is visible as a U-shaped negative bifurcation 
(R c 7 ). Time-averaged reverse flow between the saddle point (Sfj) and the fin is evidence of recir¬ 
culation produced by the first horseshoe vortex (Fj 1 ); accompanying this vortex is the corner eddy 
(F e 1 ) which runs along the root of the fin. 

In Fig. 19, the line of low mean shear is indicated by the shaded region between the half saddles 
(S e 2 and S' e2 ) and is topologically consistent with the recirculation eddy (F e 2) which accompanies 
the second horseshoe vortex (Fp). As shown in the surface-oil flow pattern (Fig. 16), this line of 
low mean shear merges with the line of negative bifurcation (B~) to form the legs of the U-shaped 
separation. 

Figure 20 shows that, in comparison with the first and second vortices (Fj\ and Fp), the third 
vortex (Fp) and its accompanying eddy (Fp) are much smaller, and given their close proximity 
to the primary stagnation point (Sfj) they are largely embedded in the leading-edge separation 
(B\7). Since all three vortices (Fj 1, Fp and Fp) have the same rotation sense, the legs on each 
side of these horseshoe vortices eventually merge with one another by their induced velocity; the 
combined vortex legs run adjacent to the U-shaped negative bifurcation (R ( 7 ). 
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(a) 

Fj i + F e i = Enodes = 2 


JSp + ±s' n + 4s,i + ^S' el = Esaddles = 2 

(b) 

Fj\ + F e i + Fj 2 + F e 2 = Enodes = 4 


is n + \S'j ! + 45 e i + 45'j + Sj 2 + $S e2 + = E saddles = 4 

(c) 

Fj] + F e] + F /2 + F el + F /3 + F e3 = Enodes = 6 


+ ±S'j + Sj 2 + ±Se 2 + + S,- 3 + ±S e3 + = Saddles = 6 


Figure 15: Junction flow at the leading edge of the fin with side-view regional topology 
£ nodes — Esaddles = 0, after [11]. The mean flow direction is from left to right. 
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(b) Primary stagnation point U-shape separation 



Figure 16: (a) Time-averaged surface-oil flow pattern around a NACAfin [11]. (b) Streakline 
interpretation. The mean flow direction is from left to right. For details of the flow 
shown in the dashed box, see Figs. 17-20. 
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Figure 17: An enlarged view of Fig. 16 showing the time-averaged surface-oil flow pattern at the 
junction of a NACA fin [11 ]. 
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Section A-A 



Section A-A 

See Fig. 15(a) for details 

Plan view 

N e i = £ nodes = 1 

Sfj = £ saddles = 1 


Figure 18: Streakline interpretation of Fig. 17 showing plan-view 
£ nodes — £ saddles = 0 with one horseshoe vortex. 
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Section A-A 



Section A-A 

See Fig. 15(b) for details 

Plan view 

N e i + N e 2 = £nodes = 2 

Sfj + .S', 2 = £ saddles = 2 


Figure 19: Streakline interpretation of Fig. 17 showing plan-view regional topology 
£ nodes — £ saddles = 0 with two horseshoe vortices. 
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Section A-A 



Section A-A 

See Fig. 15(c) for details 

Plan view 

N e i + N e 2 + N e -i = £nodes = 3 

S fj + S' e 2 + S' e3 = £ saddles = 3 


Figure 20: Streakline interpretation of Fig. 1 7 showing plan-view regional topology 
£ nodes — £ saddles = 0 with three horseshoe vortices. 
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7 Flow at the Tip of the Fin 

Flow topology of the tin is inferred from available images of surface-oil flow visualisation and 
smoke visualisation for similar NACA profiles at yaw (\|f ~ 10°-12°) [13, 14]. In Figs. 21 and 22, 
the oil streaks track the time-averaged surface-flow bifurcations. In Fig. 23(a), the smoke filaments 
(illuminated by a 2 mm thick laser-light sheet) capture the instantaneous roll-up. In the simplest 
interpretation given in Fig. 23(b), the tip vortex (F/ t ) rolls up on the leeward side; the adjacent 
half-saddles (Sf t and S'p) satisfy Eq. (3) and are topologically consistent with the positive- and 
negative-bifurcation lines (Bp and BJ t ) observed on the fin (Fig. 22). As the size of the tip vortex 
(Ff t ) increases in the chord-wise (x) direction, the distance between the positive- and negative- 
bifurcation lines (Bp and BJ t ) also increases. The departure of the tip vortex from the trailing 
edge of the fin is accompanied by a slight upward deflection of the positive-bifurcation (Bp) line. 

A closer inspection of Fig. 23(a) reveals a pair of counter-rotating eddies (F te and F( e ) near the 
root of the tip vortex (. Ff t ). These eddies are much smaller than the tip vortex and they are rapidly 
entrained into the tip vortex [14]. In the time-averaged surface-flow pattern, the counter-rotating 
eddies are embedded in the larger “footprint” of the tip vortex, i.e. the Bp line — the region of 
low mean shear indicated by the shaded area in Fig. 22. In the alternate inteipretation given in 
Fig. 23(c), the topology of the fin-tip vortex includes the pair of counter-rotating eddies (F te and 
F’ e ). In the region of low mean shear indicated by the shaded area (Fig. 23(c)), the accompany¬ 
ing saddle points (S'p, S te and S' te ) appear as very fine alternating local positive- and negative- 
bifurcation lines. These alternating bifurcation lines converge as the tip vortex separates from the 
fin into the downstream flow. 

Figure 24 shows that, at the trailing edge of the fin, a change in the yaw angle and/or the (chord) 
Reynolds number can produce more or less intense mixing of the smoke filaments, but this does 
not alter the basic topology, i.e. the cross-stream would separate at the tip and roll up to produce a 
stable focus [13, 14]. From kinematic principles (e.g. [3]), it is possible to deduce some features 
of the tip vortex downstream of the fin. For example, as the vortex/focus separates from the fin, 
the two half-saddles (i.e. Sp and S'p) connect to form a full saddle as shown in Fig. 25. Also, as 
the focus (Fft) interacts with the surrounding fluid, the vortex stretches and increases in size due 
to the effect of viscosity. A stretched vortex may contain “limit cycles”; Fig. 25 shows an example 
of a limit cycle which represents a change from inflow (stable focus) to outflow (unstable focus); 
the unstable focus, once established, extends downstream to infinity. 


8 Vortex Skeleton and Circulation 


Earlier sections in this report have identified various features of a submarine flow. These include 
the hull vortices, the fin-junction vortices and the fin-tip vortex. From this, it is possible to con¬ 
struct a vortex skeleton of the submarine flow by using the first vortex law of Helmholtz [15, 16]. 
The law states that “the circulation around a vortex tube 4 is constant ” [16]. The inference from 
this is that the vortex tube (or line) must either form a closed loop, extend to infinity or terminate 

4 A vortex tube is a set of vortex lines passing through a simply-connected surface in space, a vortex line is a line 
which is everywhere tangent to the local vorticity vector [16]. 
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Leeward 



Dimensionless chord length 


Figure 21: Time-averaged surface-oil flow pattern on a NACAfin [13] (in 3rd-angle orthogonal 
projection): \|/ ~ 10°-12°. The mean flow direction is from left to right. 
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0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


Dimensionless chord length 


Figure 22: Streakline interpretation of Fig. 21 showing top-view and side-view regional topology 
£nodes — Esaddles = 0. Shaded area indicates region of low mean shear. 
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(a) 

Cross-section smoke visualisation 

(b) 

Ff t = £ nodes = 1 

ISp + \S' ft = Esaddles = 1 

(c) 

Ff t + F le + F' le = E nodes = 3 

kSft + {S' ft + \S te + fS' te + S" = Esaddles = 3 


Figure 23: Development of a NACA fin-tip vortex for \|/ ~ 10°-12°. (a) Smoke filaments in air 
[14]. (b) Interpretation of regional topology (X nodes — Esaddles = 0) with the aid of 
surface streaklines from Fig. 22. (c) Alternate interpretation with a pair of counter¬ 
rotating eddies (F te and F' te ) at the root of the fin-tip vortex (Ff t ). 
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Figure 24: Topology interpretation of the tip vortex aft of a NACA fin; smoke filaments in air [14] 
with (a) \\t = 4° at the chord Reynolds number Re c = cyLL/v = 3000, (b) \|/ = 12° at 
Re, = 3000 and (c) \\t = 12° at Re c = 7000. The mean flow direction is out of the page. 


Top view Helical streamline Vortex tube with increasing 



Figure 25: Kinematics of the separated fin-tip vortex (Ff t ) in the very near field. In cross-section, 
the flow is viewed from downstream of the fin. 
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at solid surfaces 5 . For simplicity, the following analysis excludes vortices which may be produced 
by the control surfaces and/or support structures for the submarine model. 


Figure 26 provides an interpretation of the vortex skeleton of the separated flow for small angles of 
yaw (0 ° < \|/ < 1 8 °). The cross-sections of the vortex skeleton, interpreted from standard topology 
rules, are shown in Fig. 27. On the leeward side of the submarine, vortex lines form in the vicinity 
of the nose to produce a pair of counter-rotating hull vortices. For the upper-hull and lower-hull 
vortices, their centroids are at the vertical locations z+ hc and z c hb , respectively. At the fin, the 
leeward flow may be interpreted as a supeiposition of the free-stream and a U-shaped vortex 6 . 
This U-shaped vortex (line) consists of three segments: (i) a bound vortex which spans the height 
of the fin about the aerodynamic center 7 of the fin, connected to (ii) a vortex leg at location z~j t 
along the tip of the fin and (iii) a vortex leg at location z+ hc along the casing. 


Figures 26 and 27 show that, on the casing (at location zt h ), the rotation sense of the leg of the 
U-shaped vortex produced by the fin is opposite to that of the upper-hull vortex stretched from the 
nose, and so there is some cancellation in circulation. Note that fin-junction vortices are present 
on the casing (at starboard z~j s and port f p ) but they are much smaller than the hull vortices 
and the fin-tip vortex. For the fin-junction vortices, their size would scale approximately with the 
displacement thickness 8 of the boundary layer on the casing, i.e. ~ L/\jRei_. 

Dimensional analysis suggests that circulation (T in units of length 2 /time) of the flow around the 
submarine at yaw can be expressed as the dimensionless parameter: 


K = 


r m LL sin(\|/) 


for sin(\|/) / 0, 


(4) 


where r m is the maximum radius of the bare hull and LL sin(\|/) is the cross-stream velocity. Aft of 
the hull (or at the tail plane), it is assumed that all vorticity is shed from the submarine and gets 
wrapped up or concentrated into these identifiable longitudinal vortices in the flow field: 


the hull vortices ^ 
the fin-tip vortex j 
and the fin-junction vortices 


casing 

base 


K + = uic 

hc ^.OoSinty)’ 
K - = Ut 

hb r m il sin(\|/) ’ 


kT = 


f* r m lL sin(\|/)’ 


rt 


I starboard xt = — ‘-A —-, 
fi r m UL, sm(\|/) 


1 


port 


IN/V. — 


fp r,„tLsin(\|/)’ 


(5) 

( 6 ) 
(V) 


where K+>0 and K <0. In cross-section (Fig. 27), positive circulation is obtained by integrating 
over area of positive vorticity 9 , T + (x) = ff A oo/ cl A > 0, and negative circulation is obtained by 
integrating over area of negative vorticity, r~(x) =/’/’, (O f cl A < 0. 

5 Section 4, for example, has shown that a sink focus (F/ — a termination point of a vortex line) on the fin is 
topologically possible. 

6 Surface-flow visualisation has shown evidence of a U-shaped vortex on the leeward side of the fin, see Section 4. 

7 The aerodynamic center point of a symmetrical NACA aerofoil is approximately 25% of the chord from the leading 
edge. For an angle of attack up to 18 °, the moment coefficient about the quarter-chord point is approximately zero [5], 

8 For example, the displacement thickness of a boundary-layer flow past a flat plate parallel to the flow is given by 
8* = 1.721 x/\/Re x [17]; the boundary layer is thin provided that the Reynolds number is large. 

9 The body-axis (x) component of vorticity from the velocity field is defined as cq v = —J/. Vorticity may be 

calculated based on Stokes theorem, cq v = T/A = (1/A) § c U yz ■ dl, where the circulation T around the boundary C of an 
area A is equal to the integral of vorticity over the area A. 
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Figure 26: Schematic diagram of vortex lines and circulation around the submarine. The 
mean flow direction is from left to right. The submarine yaw angle is in the range 
0° < \J/ <18°. The maximum radius of the bare hull is defined as r m . 
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Figure 27: Schematic diagram of cross-stream around the submarine. The flow is viewed from the 
downstream end of the hull; K + >0 and K - <0 denote counter-clockwise and clockwise 
rotations, respectively; zt> 0 and z c < 0 denote the vertical locations of the centroid 
of vortices in the upper hull and the lower hull, respectively. 


An inteipretation of circulation around the submarine can be obtained by using the second vortex 
law of Helmholtz (also known as Kelvin’s theorem) [15, 16]. The law states that ‘‘the circula¬ 
tion around a material loop is time-independent, provided the fluid is inviscid, only subject to 
potential body forces, and its pressure is a function of density alone” [16]. The inference from 
this is that an initially irrotational flow (given the three constraints) generates no net circulation; 
if £r(x, y, z, t= 0) = 0, then £T = 0 for all time (?) throughout space. If a submarine is accel¬ 
erated from rest in an initially irrotational flow, the net circulation of the flow therefore would 
remain zero: 

ir(r) =0 =>- £k(x) = K+ (x) + K~ (x) + K^(x) + K+^x) + Kj p (x) + £K imag (x) = 0, (8) 

where £K imag (x) takes into account the circulation of all mirrored images inside the submarine 
geometry. Note that the fin-junction vortices on the casing are part of the same horseshoe vortex, 
and so K /v (x) + K^(x) = 0. At the tail plane x = L, Eq. (8) may be simplified to 

Lk(L) = K+ (L) + K hb (L) + K+(L) = 0 (9) 

since £ K miag (L) = 0. In Eq. (9), kT. is the residual circulation due to interactions between the 
fin U-shaped vortex leg on the casing and the upper-hull vortex stretched from the nose, where 
K /(( < Kf t . In the absence of the fin, there is no U-shaped vortex, and so K = 0 and Eq. (9) would 
reduce to KjJ" c = — Kf b . A schematic diagram of the distribution of circulation around the submarine 
is given in Fig. 26. 

To summarise, the constructed vortex skeleton of the flow around the submarine at yaw consists 
of longitudinal vortices which are stretched downstream to infinity. On the upper hull, a bound 
vortex connects the longitudinal vortices on the leeward side of the fin (Fig. 26). If the submarine 
is accelerated from rest in an initially irrotational flow (with the constraints indicated by Kelvin’s 
theorem), the net circulation of these vortices around the submarine would remain zero for all time 
throughout space. 
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9 Sway-Force and Yaw-Moment Equations 

An understanding of the flow structure (or vortex skeleton) and circulation around the submarine 
is useful because this can assist with the interpretation of the hydrodynamic forces. In this section, 
the sway forces are calculated on a part-by-part basis. On the hull, the sway force due to the 
longitudinal vortices is determined by using the “hydrodynamic-impulse” method described by 
Jeans et al. [18, 19, 20]. On the fin, the sway force due to the bound vortex is obtained separately 
by applying the first vortex law of Helmholtz [15, 16]. Where these sway forces act about the 
mid-ship (or the center of buoyancy) of the submarine, this gives an estimate of the yaw moment. 


9.1 Load on the Hull Form 

Jeans et al. [18, 19, 20] have shown that, from the slender-body theory [21, 22], the force distribu¬ 
tion parallel to the cross-stream (y) direction of a hull at yaw is 

Fy,/, = pLLcos(\|/) f di yh {x ) (10) 

Jo 

where i y j,(x) is the hydrodynamic impulse in the cross-stream (y) direction of the hull with cross- 
section area A/,(x): 

/ OO p oo 

/ (cOjcz) dydz - LL sin(\|/) A h (x ). (11) 

oo J —oo 

On the right side of Eq. (11), the first term is a moment of vorticity. Aft of the hull, the second 
term is zero given that A/,(L) = 0. 

Since the moment of vorticity depends on the vertical distribution of longitudinal vortices (i.e. 
tOyz), it is possible to model Eq. (11) by using the vortex skeleton described in Section 8. To begin, 
let the hydrodynamic impulse be defined by the system of vortices at the tail plane x = L: 

= r m LL > sin(\|/) x E (k Zc,aft) 

= r m LL sin(\|/) x ^K^.z^ c aft + K hb z c to aft + K /^/f,aft) ^ 0, (12) 

where the circulation parameter K and the vertical centroid location z c ,ait of each tail-plane vortex 
are assumed to be roughly constant for limited yaw angles in the range 0° <\|/< 18°. Eq. (12) 
excludes the much smaller fin-junction vortices because they are part of the same horseshoe 
vortex formed on the hull casing and their moments of vorticity would approximately cancel 
(i.e. K+z+ fi + K fpZ^f p — 0, see Figs. 26 and 27). For the hull vortices and the fin-tip vortex, they 
are asymmetrically distributed about the hull and their circulations are different, and so the net 
moment of vorticity due to these vortices would not equal to zero. If the net moment of vorticity is 
zero, there would be no sway force acting on the submarine. 

From integrating Eq. (10) and substituting Eq. (12) into the expression, this yields 

Fy,/, = pUir m cos(t|/)sin(t|/) X + ^c,hb,aft + K+ z+ /f aft ) . (13) 

The mid-ship yaw moment due to this force may be expressed as 

M, a ,„ — Kj t x YyjJ . m , (14) 

where L m is the distance from the tail plane to mid-ship and K] x is a coefficient of proportionality. 
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9.2 Load on the Fin 

In Section 8, the wake of the fin is interpreted as a superposition of the free-stream and a U-shaped 
vortex. In reality, the wake of the fin would contain infinitely many instantaneous U-shaped- 
like vortex lines. These vortex lines are spread along the fin to form a vertical vortex sheet with 
local strength F/i,(z ); the sway force acting on the aerodynamic center of the tin is given by the 
expression: 

rh c +hf 

F y j = plL I F fb {z)dz, (15) 

J h c 

where h c is the height of the casing and hj is the height of the fin measured from the casing 
(Fig. 26). By applying the first vortex law of Helmholtz [15, 16], the circulation produced by the 
fin bound vortex is equal to the circulation produced by the fin-tip vortex, i.e. 


F fb = T/O 


(16) 


and so Eq. (15) may be recast as 


Fyj = pLL 2 r m sin(\|/) x K f K^ t h f , (17) 

where Kf is a coefficient of proportionality and k| ; is defined by Eq. (6). The mid-ship yaw 
moment due to this force may be expressed as 

f m F y f X [L L m ^C/a) 5 (1^) 

where L c / 4 = Lf + Cf/4 is the distance from the submarine nose to the quarter-chord point of the 
fin (see Fig. 26); for the Joubert geometry, Lf/L = 0.31, Cf/L = 0.16 and L c / A /L = 0.35. 


9.3 Load on the Overall Geometry 


The total sway force acting on the submarine is obtained by adding the sway forces given by 
Eqs. (13) and (17): 


Fy — F^/, + F yj 

= p U^L r m x [K\ cos(\|/) + K f K 2 ] sin(\|/) (19) 

and the total yaw moment about the mid-ship is obtained by adding the yaw moments given by 
Eqs. (14) and (18): 


pci L 2 /•„ 


X 


'"Mcos( ¥ )+ll- l 


'c/A 


K f K 2 


sin(\|/). 


( 20 ) 


where the non-dimensionalised moments of vorticity are given by the expressions: 

Ki = 


he 


7 + 7 7+ 

he,ait , — Jib,aft , + ^c, ft,aft 

—+%—+4—. 


, he 

*2 = 4f- 


( 21 ) 

( 22 ) 
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From non-dimensionalising Eqs. (19) and (20) by the free-stream velocity LL and the submarine 
length L, this gives the sway-force and mid-ship yaw-moment coefficients: 


1 


, 2 t 2 


C Fy = Fy [ -pUjL 


K\ , K,K 2 
— cos(v) + 


sin(\|/), 


CM ; — 7vl ; 

C/;j K/jKl 

~L R 


1 


2 r 3 


2 P U ~ L 


cos( ¥ )+ 


74 \KfK 2 


LJR 


sin(\|/), 


(23) 


(24) 


where R=L/(2r m ) is the slenderness ratio of the bare hull; at mid-ship, L m /L = 1/2. 


10 Circulation, Force and Moment: A Review of Data 

This section is set out to find the model coefficients (K\, K 2 , Kf and K/, in Eqs. (21)-(24)) for the 
Joubert submarine geometry. To determine the model coefficients, this requires force and moment 
data (Cp v , Cm.) and circulation data (k, z.<). 

At the time of writing, the author knows of two cases where circulation, force and moment data 
have been published; they are from RANS CFD 10 simulations for bare hulls “DRDC-STR” and 
“Series 58-4621” [18, 19, 20]. On the basis of geometric similarity (i.e. hulls with similar slen¬ 
derness ratios, see Fig. 28), it seems likely that these bare hulls and the Joubert hull would also 
produce similar circulation of cross-stream vortices. Given that the Joubert fin has the shape of 
a rounded-tip NACA-0015 aerofoil, circulation of the fin-tip vortex may be estimated by using 
available data for a stand alone NACA-0015 aerofoil 11 . 


10.1 Scaling of Circulation 

Figures 29 and 30 provide a review summary of circulation data for the bare hulls DRDC-STR 
and Series 58-4621 and for the aerofoil NACA-0015. 

The bare-hull parameters: 


r 

r m LL sin(\|/) 


2 R 



tan(\|/) x r* Ap , 


Zc 

L 



tan(t|/) x z* Ap 


(25) 

(26) 


10 The RANS CFD (Reynolds averaged Navier Stokes computational fluid dynamics) are calculated using the SST 
(shear stress transport) “k-co” turbulence model, see [18, 19, 20]. 

11 Note that, for a fin-and-hull combination, the sway force may differ from that of a stand alone fin because of 
interference effects, e.g. [25, 26]. For ease of analysis, possible fin-and-hull interference effects are absorbed by the 
coefficient of proportionality Kf introduced by Eq. (17). 
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Hull form 

R = L/(2r m ) 
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Figure 28: Examples of generic hull forms: (a) Defense Research and Development Canada — 
Static Test Rig, (b) Naval Surface Warfare Center Carderock Division Series 58-4621, 
and (c) Defence Science and Technology Group — the submarine shape by Joubert. 
Geometric parameters are defined in Fig. 26; n/a = not applicable. 


UNCLASSIFIED 


35 


















DST-Group-TR-3177 


UNCLASSIFIED 


are plotted as functions of body length (x/L) in Fig. 29 and as functions of yaw angle (\\i) in 
Fig. 30, where the values for the dimensionless groups: 

r 


n = 


V,Ap 


Uoo sin(\|/) ^Aptan^) ’ 

Zc 

^A p tan(\\i) 


(27) 

(28) 


are taken from RANS CFD simulation [20]. For both DRDC-STR and Series 58-4621, the non- 
dimensionalised planform area, A p /L 2 , is approximately 0.1. 


In Figs. 29 and 30, the tin parameters are defined as 

r _ c f 


K = 


r m IJL sin(\|/) r m sin(\|/) 


x r* 


A Cf 

— = T X Zc C’ 

L L 

where Cf is the chord length; the values for the dimensionless groups: 

r* = r 

C UooCf’ 

* _ Zc 

'■c.c 

Cf 


(29) 

(30) 


(31) 

(32) 


are taken from various wind-tunnel experiments [7, 27, 28] and RANS CFD simulation [29]. In 
Figs. 29 and 30, the NACA-0015 aerofoil data (k, z c /L) are plotted using the ratios Cf/L = 0.16 and 
Cf/r m = (cf/L) x2R = 2.34 for the Joubert submarine geometry (see Fig. 28). The leading edge of 
the Joubert fin is at x/L = 0.31. 


10.2 Interpretation of Circulation Data 


In Figs. 29 and 30, K for the bare hulls is the circulation magnitude of the counter-rotating vortices 
(i.e. K/; C = K hb)\ for the fin, K is the total strength of the tip vortex (i.e. Kj,). For each structure, its 
vertical centroid location z c is measured from the body (x) axis of the hull (see Figs. 26 and 27). 

Figure 29(a) shows that circulation increases with increasing streamwise distance from the lead¬ 
ing edge of both the hull and the fin. The rise in circulation is most likely due to the summation 
of circulation contributed by an increasing number of vortex lines formed along the submarine. 
Figure 31 shows a schematic diagram of possible distribution of vortex lines, where (finer) up¬ 
stream vortex lines eventually roll up or grouped together to form (thicker) downstream vortex 
tubes. 


From inspection of Fig. 29(b), it is possible to infer the location of the vortices at the tail plane. 
For example, the vertical distance between the centroids of the hull vortices does not exceed the 
maximum diameter of the hull, i.e. 


*■<?, he, aft <- c, hb, aft 


= lh 


1 


0 < Y/, < 1, 


L " R’ 

and the vertical location of the tip vortex does not exceed the height of the fin, i.e. 

1 


■7+ 

aft 


hf 

“ Y/ l 2 R + T 


0 < Y/ < 1, 


(33) 


(34) 
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(a) 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Body length, x/L 


Figure 29: Circulation and location of vortices along the submarine. Eqs. (25)-(32) are used to 
plot the data; bare hull [20]; fin “ -” [7], “A” [27]. 
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Figure 30: Circulation and location of vortices as functions of \|/. Eqs. ( 25)-(32 ) are used to plot 
the data; bare hull [20]; fin “•”[?], “A" [27], “sy” [28], “U” [29]. 
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Bound vortex 


Surface 


Fin-tip 

vortex 
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fine-scale 'eddies 


Hull vortices F hl 



Figure 31: Possible distribution of vortex lines and circulation around the submarine, after 
Fig. 26. The mean flow direction is from left to right. For clarity, the junction vor¬ 
tices around the fin are not included. 
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where y* and y f are coefficients of proportionality for the hull vortices and the tin vortex, respec¬ 
tively. For a hull with an approximately round cross-section, this gives z+ hc aft — —z~ hh aft , and so 
by substituting Eqs. (33) and (34) into Eq. (21) this yields 

K 'A < 35 ) 

which relates coefficient K\ to the shape of the submarine. 

It should be noted that Eq. (35) is only an approximate for a narrow range of yaw angles. For 
example, Fig. 30 shows that with increasing yaw angle (\|t), there is a slight increase in circulation 
(k) for the vortices shed by the hull and the tin; the vertical locations of these vortices (z, ) are 
approximately constant. As yet, the present analytical model does not take into account detailed 
variations in the flow (circulation) with changing yaw angle 12 . 

10.3 Calibrating the Force and Moment Equations 

Section 9 has established that the equations for the sway-force and yaw-moment coefficients are 
of the forms: 


where 


Cf, = 

[A cos(\|f) + B] sin(\|/), 

(36) 

Cm, = 

[C cos(\|/) + D] sin(\|/), 

(37) 

A = 

R ’ 

(38) 

B = 

K f K 2 

R ’ 

(39) 

C = 

t k " a - 

(40) 

D = 


(41) 


For yaw moment about the mid-ship, L m /L = 1/2. In Fig. 28, the hull forms are shown with and 
without a fin. In the absence of a fin, Kf and Ki (and hence B and D ) are equal to zero, and K\ and 
K/, (and hence A and C) are solved simultaneously by least-squares fitting to force and moment 
(Cf v , Cm.) data. With a fin, K\ and Ki are obtained by solving Eqs. (35) and (22) respectively, and 
then the remaining coefficients Kf and K \ are solved simultaneously by curve fitting to force and 
moment (Cf v , Cm.) data. 


10.3.1 The “DRDC-STR” and “Series 58-4621” Bare Hulls 

In Figs. 32 and 33, the force and moment coefficients are plotted as functions of yaw angle (\|r) for 
the bare hulls DRDC-STR and Series 58-4621. The data are from measurements in a wind tunnel 

l2 At very small yaw angles (0° < V|/^C 10°), the flow is nearly symmetrical about the mirror plane of the submarine 
with negligible leeward vortices produced along the hull (see Section 4). At larger yaw angles (\|/>10°), strong circu¬ 
lation produces a reattachment node on the leeward face of the fin; further increasing the yaw angle can produce greater 
separation in the vicinity of the nose leading to the Werle-Legendre horn vortices (see Section 5). 
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and in a towing tank, and from RANS CFD simulations, see [20]. Inspection shows that, due to 
scatter, there is little discernible difference between the data for DRDC-STR and Series 58-4621; 
the data (<, ◄, o, ►) in the range 5°<\|/<15° fall on 

Cp Y = 0.015 cos(\|/) sin(\|f), (42) 

Cm- = 0.011 cos(\|/) sin(\|/), (43) 


with root-mean-square (r.m.s.) errors of 26% and 8%, respectively. 

From Eq. (9), £k = K /jr + k J lb = 0 since Kj t = 0 in the absence of the fin. So, by substituting 
Eq. (38) into Eq. (35) with = 0 and rearranging the expression, this gives the circulation of the 
bare-hull vortices: 


K hc~ K hb _ AxR 2 

2 y» 



~ K h h = 


Ax R 2 
Th 


(44) 


where A = 0.015 from Eq. (42). 


Inspection of Fig. 30(b) shows that, at x = L, the vertical location of the bare-hull vortices is 
approximately one-third of the maximum radius of the bare hull, i.e. y ), = 0.34 for DRDC-STR 
and y/ ; = 0.28 for Series 58-4621. Solving Eq. (44) yields = 0.015 x8.75 2 /0.34 = 3.4 

for DRDC-STR and K^. = —k J lb = 0.015 x7.34 2 /0.28 = 2.9 for Series 58-4621. These predicted 
values of K are in agreement with the bare-hull simulation data (◄, ►) shown in Fig. 30(a) for the 
range 5 ° < \|/ < 15°. 


10.3.2 The “Joubert” Hull Form without a Fin 

Figures 34 and 35 show the force and moment coefficients plotted as functions of yaw angle (\|/) 
for the Joubert submarine 13 starting with a bare hull, then adding the casing and the fin. In the 
range 5 ° < t|/ < 15 °, the Joubert bare-hull data (o) fall on 

Cf y = 0.036 cos(\|/) sin(\|/), (45) 

Cm. = 0.014 cos(\|/) sin(\|/), (46) 

with r.m.s. errors of 20% and 7%, respectively. With the casing, the force and moment coefficients 
are slightly larger and the experimental data (+) fall on 

Cf y = 0.042 cos(\|/) sin(\|/), (47) 

Cm. = 0.017 cos(\|/) sin(\|/), (48) 

with r.m.s. errors of 21% and 7%, respectively. Note that the casing adds to the height of the 
hull and so it seems likely that the vertical distance between the hull vortices, i.e. y ), as defined in 
Eqs. (33) and (44), would also increase with the height of the hull. A plot of Eq. (44) in Fig. 36(a) 
shows the effect of increasing y/, on the circulation (k f, K /ih ) for the Joubert bare hull without the 
casing (A = 0.036) and with the casing (A = 0.042). 

13 The measurements are for Rei = 5.2x 10 6 tested in the Defence Science and Technology Group low-speed wind 
tunnel [23, 24]. They are obtained by using a 6-component strain-gauge balance fitted inside the submarine via a 
single-pylon support on a turntable. The measurement uncertainties for Cf v and Cm. are 0.15xl0~ 3 and 0.03 xl0~ 3 , 
respectively [24]. 
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Figure 32: Force and moment coefficients as functions of yaw angle for generic bare hulls; a 
cur\’e fit ofEqs. (23) and (24) on the review data from [20], 
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Figure 33: Logarithmic-linear plot of Fig. 32. 
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Figure 34: Force and moment coefficients as functions of yaw angle for the generic Joubert hull 
form; a curve fit ofEqs. (23) and (24) on the experimental data from [23, 24]. 
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Figure 35: Logarithmic-linear plot of Fig. 34. 
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Figure 36: Summary of least-squares fitting to force and moment data for the Joubert hull form 
using Eqs. (22), (35) and (36)-(41). (a) Bare hull with and without casing, (b, c) Bare 
hull with casing and fin. R = 1 .3, L c u= 0.35, hf/L = 0.08 and L m /L= 1/2. 
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10.3.3 The “Joubert” Hull Form with a Fin 

To curve fit the fin-appended data (o) in Figs. 34 and 35, the process begins by solving K 2 in 
Eq. (22) and K\ in Eq. (35). This requires information on the fin-tip vortex (k/, Jj) and the hull 
vortices K^, and y/,); the information is obtained as follows. 

1. Inspection of Fig. 30 shows that, for the range 5°< \|/ < 15°, the fin-tip vortex has ap¬ 
proximately constant circulation k/ ; = 3 and is located at about three-quarter of the dis¬ 
tance from the center-line of the hull to the top of the fin, i.e. Jf = 3/4. By using Eq. (22), 
K 2 = k+ x (hf/L) = 3 x 0.08 = 0.240. 

2. Note that no alteration is made to the lower hull, therefore circulation on the lower hull is 
assumed to be unchanged, and so Kj ih = —A x R 2 /ji, from Eq. (44), where A = 0.042 from 
Eq. (47). 

3. Since the flow field on the upper hull is dominated by the fin-tip vortex, much of the circu¬ 
lation on the upper hull is contributed by this vortex with the constraints: 

/tv/? 2 /tv/? 2 

Kp > <C > 0 and -^-q—<Y/ ! <——. (49) 

ZK ft K ft 

where A = 0.042, R = 7.3 and k^ = 3 specify the range 0.373 < y/, < 0.746. 

4. From Eq. (9), the residual circulation on the upper hull is determined as a function of y h, 
that is, K/+ = -K hh - Kj t = (0.042 x R 2 /y h ) - 3, see Fig. 36(b). 

5. By using Eq. (35), K\ is solved as a function of ysee Fig. 36(c). 

6. Once K\ and K 2 are known from the above steps, the remaining coefficients Kj and K\, as 
functions of y ) x are obtained by least-squares fitting Eqs. (36) and (37) to the force-and- 
moment data, see Fig. 36(c). 

7. From Eqs. (38)-(41), the model coefficients (A, B, C and D) are determined as functions of 
y/„ see Fig. 36(c). 


To summarise with an example say y> 7 = 1/2, this gives a lower-hull circulation — K /(/) =A x R 2 /jh = 
0.042 x 7.3 2 /( 1/2) =4.5. A fin-tip vortex of strength = 3 leaves on the upper-hull casing a 
residual circulation K hc = —k J ib — Kj t = 4.5 — 3 = 1.5. Hence, substituting the values K J ib = —4.5, 
Kyy = 3, K / j c = 1.5, y i, = 1/2 and yy = 3/4 into Eqs. (35) and (22) yields K\ =0.540 and Ki =0.240. 
By least-squares fitting the data (o, Figs. 34 and 35) in the range 5 ° < V|/ < 15 °, 

C Fy = [0.074 cos(\|/) + 0.0192] sin(\|/), (50) 

Cm, = [0.020 cos(\|/) + 0.0029] sin(\|/), (51) 

withr.m.s. errors of 9% and 3% respectively, the remaining coefficients C/ =0.584 and K/, = 0.546 
are solved. 
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11 Concluding Remarks 

A model of the flow around a generic conventional “Joubert” submarine hull form at pure yaw has 
been constructed from surface-streamer visualisation and flow-topology concepts. The purpose of 
this work is to gain some understanding of the large-scale (vortical) structures produced by the 
hull which can affect the motion control, manoeuvring and signature of the submarine. 

The flow-topology model (e.g. Fig. 31) shows the presence of a pair of counter-rotating longitu¬ 
dinal vortices on the leeward side of the hull. The inclusion of a fin (or sail) on the upper hull 
further produces a U-shaped vortex. This U-shaped vortex consists of three segments: (i) a bound 
vortex which spans the height of the fin, connected to (ii) a vortex leg along the tip of the fin and 
(iii) a vortex leg along the upper hull. Just above the upper hull, the rotation sense of the vortex leg 
arising from the fin is opposite to that of the vortex leg arising from the nose, and so these vortex 
legs act to cancel each other (the resulting structure F), c is shown in Fig. 31). 

An analytical treatment of the topology model has provided an initial step towards developing 
equations which relate the sway-force and yaw-moment coefficients to the hull-form geometry 
and the circulation of the surrounding flow. The model equations (36) and (37) take into account 
the fin location and the slenderness ratio of the hull. They can be calibrated by using circulation, 
force and moment data from experiments and/or computational-fluid-dynamics (CFD) modelling, 
and therefore can serve to cross-check/validate experimental and CFD data. 

Future refinements may take into account other factors such as casing geometry, control surfaces, 
experimental arrangement of sting or pylon support(s), submarine propulsion and the Reynolds 
number to obtain general equations for sway force and yaw moment. 


12 Acknowledgements 

Thanks go to Chetan Kumar and Paul Jacquemin for assisting with the flow-visualisation exper¬ 
iment, to Gregory Seil and Lincoln Erm for refereeing this document, and to Howard Quick for 
providing the experimental force and moment data for the generic conventional Joubert submarine 
geometry. The financial support from the SEA 1000 project is gratefully acknowledged. 


References 

1. Joubert, P. N. (2004). Some aspects of submarine design part 1: hydrodynamics. Technical 
Report 1622, DSTO. 

2. Joubert, P. N. (2006). Some aspects of submarine design part 2: shape of a submarine 2026. 
Technical Report 1920, DSTO. 

3. Hunt, J. C. R., Abell, C. J., Peterka, J. A., and Woo, H. (1978). Kinematical studies of the flows 
around free or surface-mounted obstacles; applying topology to flow visualisation. Journal of 
Fluid Mechanics, 86:179-200. 

4. Tobak, M. and Peake, D. J. (1982). Topology of three-dimensional separated flows. Annual 
Review of Fluid Mechanics, 14:61-85. 


48 


UNCLASSIFIED 


UNCLASSIFIED 

DST-Group-TR-3177 

5. Goett, H. J. and Bullivant W. K. (1938). Tests of N.A.C.A. 0009, 0012, and 0018 airfoils in the 
full-scale tunnel. Report 647, NACA. 

6. Erm, L. R (2003). Calibration of the flow in the extended test section of the low-speed wind 
tunnel at DSTO. Technical Report 1384, DSTO. 

7. Kumar, C., Manovski, P., and Giacobello, M. (2012). Particle image velocimetry measurements 
on a generic submarine hull form. In Proceedings of the 18th Australasian Fluid Mechanics 
Conference, paper 188, Launceston, Australia. 

8. Anderson, B., Chapuis, M., Erm, L., Fureby, C., Giacobello, M., Henbest, S., Jones, D., 
Jones, M., Kumar, C., Liefvendahl, M., Manovski, P., Norrison, D., Quick, H., Snowden, A., 
Valiyff, A., Widjaja, R., and Woodyatt, B. (2012). Experimental and computational investigation 
of a generic conventional submarine hull form. In 29th Symposium on Naval Hydrodynamics, 
Gothenburg, Sweden. 

9. Fureby, C., Anderson, B., Clarke, D., Erm, L., Giacobello, M., Henbest, S., Jones, D., 
Nguyen, M., Johansson, M., Jones, M., Kumar, C., Lee, S.-K., Manovski, P., Norrison, D., 
Patterson, K., Seil, G., Woodyatt, B., and Zhu, S. (2014). Experimental and numerical study of a 
generic conventional submarine at 10-degrees yaw. In 30tli Symposium on Naval Hydrodynam¬ 
ics, Hobart, Australia. 

10. Ashok, A. and Smits, A. J. (2013). The turbulent wake of a submarine model in pitch and 
yaw. In 51th A1AA Aerospace Sciences Meeting, paper 2013-1121, Grapevine, Texas, USA. 

11. Simpson, R. L. (2001). Junction flows. Annual Review of Fluid Mechanics, 33:415^443. 

12. Ole men, S. M. and Simpson, R. L. (1995). An experimental study of a three-dimensional 
pressure-driven turbulent boundary layer. Journal of Fluid Mechanics, 290:225-262. 

13. Chow, J., Zilliac, G., and Bradshaw, P. (1997). Turbulence measurements in the near field of 
a wingtip vortex. Technical Memorandum 110418, NASA. 

14. Giuni, M. (2013). Formation and early development of wingtip vortices. Ph.D. Thesis, School 
of Engineering. University of Glasgow, United Kingdom. 

15. Helmholtz, H. (1867). LXIII. On integrals of the hydrodynamical equations, which express 
vortex-motion. Philosophical Magazine Series 4, 33:226, 485-512. 

16. Green, S. I. (1995). Fluid Vortices, Kluwer Academic Publishers. 

17. Munson, B. R., Young, D. F., and Okiishi, T. H. (2002). Fundamentals of Fluid Mechanics, 

4th edition, John Wiley & Sons, Inc. 

18. Jeans, T. L. (2007). An estimation method for the normal force distribution on slender ax- 
isymmetric bodies with tapered tails. Ph.D. Thesis, Department of Mechanical Engineering. 
University of New Brunswick, Canada. 

19. Jeans, T. L., Watt, G. D., Gerber, A. G., Holloway, A. G. L., and Baker, C. R. (2009). High- 
resolution Reynolds-averaged Navier-Stokes flow predictions over axisymmetric bodies with 
tapered tails. American Institute of Aeronautics and Astronautics Journal, 47:19-32. 

UNCLASSIFIED 

49 


UNCLASSIFIED 

DST-Group-TR-3177 

20. Jeans, T. L., Holloway, A. G. L., Watt, G. D., and Gerber, A. G. (2010). A force estimation 
method for viscous separated flow over slender axisymmetric bodies with tapered tails. Journal 
of Ship Research, 54:53-67. 

21. Lamb, H. (1932). Hydrodynamics, 6th edition, Dover Publications. 

22. Jones, R. T. (1946). Properties of low-aspect-ratio pointed wings at speeds below and above 
the speed of sound. Report 835, NACA. 

23. Quick, H., Widjaja, R., Anderson, B., Woodyatt, B., Snowden, A. D., and Lam, S. (2012). 
Phase I experimental testing of a generic submarine model in the DSTO low speed wind tunnel. 
Technical Note 1101, DSTO. 

24. Quick, H. and Woodyatt, B. (2014). Phase II experimental testing of a generic submarine 
model in the DSTO low speed wind tunnel. Technical Note 1274, DSTO. 

25. Johnson, H. S. (1953). Wind-tunnel investigation at low speed of the effect of varying the 
ratio of body diameter to wing span from 0.1 to 0.8 on the aerodynamic characteristics in pitch 
of a 45 ° sweptback-wing-body combination. Research Memorandum L53J09a, NACA. 

26. Pitts, W. C., Nielsen, J. N., and Kaattari, G. E. (1953). Lift and center of pressure of wing- 
body-tail combinations at subsonic, transonic, and supersonic speeds. Report 1307, NACA. 

27. Brich, D., Lee, T., Mokhtarian, F., and Kafyeke, F. (2004). Structure and induced drag of a tip 
vortex. Journal of Aircraft, 41:113 8-1145. 

28. Anderson, E. A. and Lawton, T. A. (2003). Correlation between vortex strength and axial 
velocity in a trailing vortex. Journal of Aircraft, 40:699-704. 

29. Seil, G. and Anderson, B. (2013). Computational Fluid Dynamics Manoeuvring Simulation 
of DSTO Generic Submarine with different Fins. Client Report 2013-0385, DSTO. 


50 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


Page classification: UNCLASSIFIED 


DEFENCE SCIENCE AND TECHNOLOGY GROUP 
DOCUMENT CONTROL DATA 


1. DLM/CAVEAT (OF DOCUMENT) 


2. TITLE 

Topology Model of the Flow around a Submarine Hull Form 


3. SECURITY CLASSIFICATION (FOR UNCLASSIFIED REPORTS 
THAT ARE LIMITED RELEASE USE (L) NEXT TO DOCUMENT 
CLASSIFICATION) 

Document (U) 

Title (U) 

Abstract (U) 


4. AUTHOR 

5. -K. Lee 


5. CORPORATE AUTHOR 

Defence Science and Technology Group 
506 Lorimer St, 

Fishermans Bend, Victoria 3207, Australia 


6a. DST GROUP NUMBER 
DST-Group-TR-3177 


6b. AR NUMBER 
AR-016-453 


6c. TYPE OF REPORT 
Technical Report 


7. DOCUMENT DATE 
December, 2015 


8. FILE NUMBER 

9. TASK NUMBER 

10. TASK SPONSOR 

11. No. OF PAGES 

12. No. OF REFS 

- 

- 

SEA 1000 


50 

29 


13. DST Group Publications Repository 

http://dspace.dstg.defence.gov.au/dspace/ 


14. RELEASE AUTHORITY 
Chief, Maritime Division 


15. SECONDARY RELEASE STATEMENT OF THIS DOCUMENT 
Approved for Public Release 

OVERSEAS ENQUIRIES OUTSIDE STATED LIMITATIONS SHOULD BE REFERRED THROUGH DOCUMENT EXCHANGE. PO BOX 1500, EDINBURGH. SOUTH AUSTRALIA 5111 


16. DELIBERATE ANNOUNCEMENT 
No Limitations 


17. CITATION IN OTHER DOCUMENTS 
No Limitations 

18. DST GROUP RESEARCH LIBRARY THESAURUS 

Submarine; Flow Visualisation; Topology; Circulation; Sway Force; Yaw Moment 


19. ABSTRACT 

A topology model constructed from surface-streamer visualisation describes the flow around a generic conventional submarine hull form 
at pure yaw angles of 0°, 10° and 18°. The model is used to develop equations for sway-force and yaw-moment coefficients which 
relate to the hull-form geometry and the flow circulation around the submarine. 


Page classification: UNCLASSIFIED 






















